Can accretion disk properties distinguish gravastars from black holes? 
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Gravastars, hypothetic astrophysical objects, consisting of a dark energy condensate surrounded 
by a strongly correlated thin shell of anisotropic matter, have been proposed as an alternative 
to the standard black hole picture of general relativity. Observationally distinguishing between 
astrophysical black holes and gravastars is a major challenge for this latter theoretical model. This 
due to the fact that in static gravastars large stability regions (of the transition layer of these 
configurations) exist that are sufficiently close to the expected position of the event horizon, so 
that it would be difficult to distinguish the exterior geometry of gravastars from an astrophysical 
black hole. However, in the context of stationary and axially symmetrical geometries, a possibility 
of distinguishing gravastars from black holes is through the comparative study of thin accretion 
disks around rotating gravastars and Kerr-type black holes, respectively. In the present paper, we 
consider accretion disks around slowly rotating gravastars, with all the metric tensor components 
estimated up to the second order in the angular velocity. Due to the differences in the exterior 
geometry, the thermodynamic and electromagnetic properties of the disks (energy flux, temperature 
distribution and equilibrium radiation spectrum) are different for these two classes of compact 
objects, consequently giving clear observational signatures. In addition to this, it is also shown that 
the conversion efficiency of the accreting mass into radiation is always smaller than the conversion 
efficiency for black holes, i.e., gravastars provide a less efficient mechanism for converting mass to 
radiation than black holes. Thus, these observational signatures provide the possibility of clearly 
distinguishing rotating gravastars from Kerr-type black holes. 

PACS numbers: 04.50.Kd, 04.70.Bw, 97.10.Gz 



I. INTRODUCTION 



The Schwarzschild solution has played a fundamental 
conceptual role in general relativity, and beyond, for in- 
stance, regarding event horizons, spacetime singularities, 
and aspects of quantum field theory in curved spacetimes. 
However, one still encounters in the literature the ex- 
istence of misconceptions, as well as a certain ambigu- 
ity inherent in the Schwarzschild solution (we refer the 
reader to jlj for a detailed review). In this context, re- 
cently a new final state of gravitational collapse has been 
proposed, denoted as gravastars (gravitational vacuum 
stars) [2J, which represent a viable alternative to black 
holes, and their properties have been extensively inves- 
tigated. These models consist of a compact object with 
an interior de Sitter condensate, governed by an equa- 
tion of state given by p = —p, matched to a shell of 
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finite thickness with an equation of state p = p. The 
latter is then matched to an exterior Schwarzschild vac- 
uum solution. The thick shell replaces both the de Sitter 
and the Schwarzschild horizons, therefore, the gravastar 
model has no singularity at the origin and no event hori- 
zon, as its rigid surface is located at a radius slightly 
greater than the Schwarzschild radius @. These config- 
urations are stable from a thermodynamic point of view. 
The issue of the dynamic stability of the transition layer 
(an infinitesimally thin shell) against spherically symmet- 
ric perturbations was considered in [3|, by constructing 
a model that shares the key features of the gravastar 
scenario. It was found that there are some physically 
reasonable equations of state for the transition layer that 
lead to stability. This latter stability analysis was further 
generalized to an anti-de Sitter or de Sitter interior and 
a Schwarzschild (anti)-de Sitter or Reissner-Nordstrom 
exterior 0]. Recently, dynamical models of prototype 
gravastars were constructed and studied [E[ . It was found 
that in some cases the models represent stable gravastars, 
while in other cases they represent "bounded excursion" 
stable gravastars, where the thin shell is oscillating be- 
tween two finite radii. In some other cases they collapse 
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until the formation of black holes. 

In addition to this, gravastar models that exhibit 
continuous pressure without the presence of infinitesi- 
mally thin shells were introduced in [|| and further an- 
alyzed in Q- By considering the usual TOV equation 
for static solutions with negative central pressure, it 
was found that gravastars cannot be perfect fluids and 
anisotropic pressures in the 'crust' of a gravastar- like 
object are unavoidable. The anisotropic TOV equation 
can then be used to bound the pressure anisotropy, and 
the transverse stresses that support a gravastar permit 
a higher compactness than the Buchdahl-Bondi bound 
for perfect-fluid stars. A wide variety of gravastar mod- 
els within the context of nonlinear electrodynamics were 
also constructed in Q. Using the F representation, spe- 
cific forms of Lagrangians were considered describing 
magnetic gravastars, which may be interpreted as self- 
gravitating magnetic monopoles with charge g. Using the 
dual P formulation of nonlinear electrodynamics, electric 
gravastar models were constructed by considering specific 
structural functions, and the characteristics and physical 
properties of the solutions were further explored. Gravas- 
tar solutions with a Born-Infcld phantom replacing the 
de Sitter interior were also analyzed in Q. 

It has also been recently proposed by Chapline that 
this new emerging picture consisting of a compact ob- 
ject resembling ordinary spacetime, in which the vacuum 
energy is much larger than the cosmological vacuum en- 
ergy, has been denoted as a 'dark energy star' [io| . In- 
deed, a generalization of the gravastar picture was con- 
sidered in by considering a matching of an interior 
solution governed by the dark energy equation of state, 
lo = p/p < —1/3, to an exterior Schwarzschild vacuum 
solution at a junction interface. Several relativistic dark 
energy stellar configurations were analyzed by imposing 
specific choices for the mass function, by assuming a con- 
stant energy density, and a monotonic decreasing energy 
density in the star's interior, respectively. The dynamical 
stability of the transition layer of these dark energy stars 
to linearized spherically symmetric radial perturbations 
about static equilibrium solutions was also considered, 
and it was found that large stability regions exist that 
are sufficiently close to where the event horizon is ex- 
pected to form. Evolving dark energy stars were explored 
in [1^ ], where a time-dependent dark energy parameter 
was considered. The general properties of a spherically 
symmetric body described through the generalized Chap- 
lygin equation of state were also extensively analyzed in 
|13|. In the context of cosmological equations of state, 
in [14j the construction of gravastars supported by a van 
der Waals equation of state was studied and their re- 
spective characteristics and physical properties were fur- 
ther analyzed. It was argued that these van der Waals 
quintessence stars may possibly originate from density 
fluctuations in the cosmological background. Note that 
the van der Waals quintessence equation of state is an 
interesting scenario that describes the late universe, and 
seems to provide a solution to the puzzle of dark energy, 



without the presence of exotic fluids or modifications of 
the Friedmann equations. 

However, observationally distinguishing between astro- 
physical black holes and gravastars is a major challenge 
for this latter theoretical model, as in static gravastars 
large stability regions exist that are sufficiently close 
to the expected position of the event horizon. The 
constraints that present-day observations of well-known 
black hole candidates place on the gravastar model were 
discussed in (l5j . The heating of neutron stars via ac- 
cretion is well documented by astronomical observations. 
A gravastar would acquire most of its mass via accre- 
tion, either as part of its birth (e.g., during core collapse 
in a supernova explosion followed by the rapid accretion 
of a fallback disk), or over an extended period of time 
after birth. Similar heating processes have not been ob- 
served in the case of black hole candidates. However, 
the absence of a detectable heating may also be consis- 
tent with a gravastar, if the heat capacity is large enough 
that it requires large amounts of heat to produce small 
changes in temperature. Nevertheless, some level of ac- 
cretion heating is unavoidable, and provides some strong 
constraints on the gravastar model. The large surface 
redshifts employed in gravastar models implies that the 
internal energy generated per unit rest mass accreted is 
very nearly c 2 , and that the radiation emitted from the 
surface of the object should be an almost perfect black 
body. The energy evolution of an accreting gravastar is 
determined by the equation dU/dt = Mc 2 — L, where 
U is the internal energy, M is the mass accretion rate, 
and L is the luminosity. Assuming that the gravastar 
is the result of a BEC-like phase transition induced by 
strong gravity the gravastar which starts at zero 
temperature and rapidly accretes a mass AmMj un will 
be heated to a temperature as observed at infinity of 
T h w 3.1 x lO 6 ™,- 2 / 3 ^ 1 / 3 (Am/m) 1/3 K, where m is the 
mass of the gravastar in solar mass units, and £ = l/lpi is 
the length scale in Planck units at which general relativ- 
ity fails to adequately describe gravity. In the case of slow 
accretion, the temperature is given by the equilibrium 

/ . \ 1/4 

value T eq w 1.7 x 10 7 [M/M Edd \ to" 1 / 4 K, where 

MEdd = 2.3 x lO _9 mM0/yr is the Eddington mass accre- 
tion rate. By comparing these theoretical predictions on 
the temperature with the upper limits of the temperature 
of the black hole candidates, obtained through observa- 
tions, one can find some limits on £. Two black hole 
candidates, known to have extraordinarily low luminosi- 
ties, namely, the supermassive black hole in the galactic 
center, Sagittarius A*, and the stellar-mass black hole, 
XTE J1118 + 480, respectively, were considered in the 
analysis [lj|. For XTE J1118+480, due to the low mass 
and low heat capacity, the value of £ is constrained to 
£ ~ 5 x 10 3 , while for Sagittarius A* values for £ in the 
range 10 4 and 10 11 are excluded. Therefore the length 
scale for the considered gravastar models must be sub- 
Planckian. Thus, if any significant fraction of the mass of 
the gravastars is due to accretion, we should see the ther- 
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mal emission associated with that phase, which is ruled 
out for Sagittarius A* and XTE J1118 + 480. 

The question of whether gravastars can be distin- 
guished from black holes at all was also considered, from 
a theoretical point of view, in [l6j], where two basic ques- 
tions analyzed were: (i) Is a gravastar stable against 
generic perturbations and, (ii) if it is stable, can an ob- 
server distinguish it from a black hole of the same mass? 
A general class of gravastars was constructed, and the 
equilibrium conditions in order to exist as solutions of 
the Einstein equations were obtained. It was found that 
gravastars are stable to axial perturbations, and that 
their quasi- normal modes differ from those of a black hole 
of the same mass. Thus, these modes can be used to dis- 
cern, beyond dispute, a gravastar from a black hole. The 
formation hysteresis effects were ignored in this study. 
In addition to this, sharp analytic bounds on the surface 
compactness 2m /r that follow from the requirement that 
the dominant e nerg y condition (DEC) holds at the shell 
were derived in [l?} • In the case of a Schwarzschild exte- 
rior, the highest surface compactness is achieved with the 
stiff shell in the limit of vanishing (dark) energy density in 
the interior. In the case of a Schwarzschild-de Sitter ex- 
terior, it was shown that gravastar configurations with a 
surface pressure and with a vanishing shell pressure (dust 
shells), are allowed by the DEC. The causality require- 
ment (sound speed not exceeding that of light) further 
restricts the space of allowed gravastar configurations. 

The crgoregion instability is known to affect very com- 
pact objects that rotate very rapidly, and that do not 
possess an horizon. A detailed analysis on the relevance 
of the ergoregion instability for the viability of gravas- 
tars was presented in [l8|, . In [l8| , it was shown that 
ultra-compact objects with high redshift at their surface 
are unstable when rapidly spinning, which strengthens 
the role of black holes as candidates for astrophysical 
observations of rapidly spinning compact objects. In 
particular, analytical and numerical results indicate that 
gravastars are unstable against scalar field perturbations. 
Their instability timescale is many orders of magnitude 
stronger than the instability timescale for ordinary stars 
with uniform density. In the large / = m approximation, 
suitable for a WKB treatment, gravitational and scalar 
perturbations have similar instability timescales. In the 
low-m regime gravitational perturbations are expected to 
have shorter instability timescales than scalar perturba- 
tions. In [19], the analysis shows that not all rotating 
gravastars are unstable, and stable models can be con- 
structed also with J/M 2 ~ 1, where J and M are the an- 
gular momentum and mass of the gravastar, respectively. 
Therefore, the existence of gravastars cannot be ruled 
out by invoking the ergoregion instability. The gravas- 
tar model was extended by introducing an electrically 
charged component in [2(j, where the Einstein-Maxwell 
field equations were solved in the asymptotically de Sit- 
ter interior, and a source of the electric field was coupled 
to the fluid energy density. Two different solutions that 
satisfy the dominant energy condition were given, and 



the equation of state, the speed of sound and the surface 
redshift were calculated for both models. The dipolar 
magnetic field configuration for gravastars was studied 
in [21[ , and solutions of Maxwell equations in the inter- 
nal background spacetime of a slowly rotating gravastar 
were obtained. The shell of the gravastar where the mag- 
netic field penetrated was modeled as a sphere consisting 
of a highly magnetized perfect fluid, with infinite con- 
ductivity. It was assumed that the dipolar magnetic field 
of the gravastar is produced by a circular current loop 
symmetrically placed at radius a at the equatorial plane. 

The mass accretion around rotating black holes was 
studied in general relativity for the first time in [22| . By 
using an equatorial approximation to the stationary and 
axisymmetric spacetime of rotating black holes, steady- 
state thin disk models were constructed, extending the 
theory of non-relativistic accretion (23|. In these mod- 
els hydrodynamical equilibrium is maintained by efficient 
cooling mechanisms via radiation transport, and the ac- 
creting matter has a Keplerian rotation. The radiation 
emitted by the disk surface was also studied under the 
assumption that black body radiation would emerge from 
the disk in thermodynamical equilibrium. The radiation 
properties of thin accretion disks were further analyzed in 
[24 , 1251 ] , where the effects of photon capture by the hole 
on the spin evolution were presented as well. In these 
works the efficiency with which black holes convert rest 
mass into outgoing radiation in the accretion process was 
also computed. More recently, the emissivity properties 
of the accretion disks were investigated for exotic cen- 
tral objects, such as wormholes (26[, and non-rotating or 
rotating quark, boson or fermion stars and brane-world 
black holes [13, H, [H, H SH [H| . The radiation power 
per unit area, the temperature of the disk and the spec- 
trum of the emitted radiation were given, and compared 
with the case of a Schwarzschild black hole of an equal 
mass. The physical properties of matter forming a thin 
accretion disk in the static and spherically symmetric 
spacetime metric of vacuum f(R) modified gravity mod- 
els were also analyzed [33] , and it was shown that partic- 
ular signatures can appear in the electromagnetic spec- 
trum, thus leading to the possibility of directly testing 
modified gravity models by using astrophysical observa- 
tions of the emission spectra from accretion disks. 

It is the purpose of the present paper to consider 
another observational possibility that may distinguish 
gravastars from black holes, namely, the study of the 
properties of the thin accretion disks around rotating 
gravastars and black holes, respectively. Thus, we con- 
sider a comparative study of the thin accretion disks 
around slowly rotating gravastars and black holes, respec- 
tively. In particular, we consider the basic physical pa- 
rameters describing the disks, such as the emitted energy 
flux, the temperature distribution on the surface of the 
disk, as well as the spectrum of the emitted equilibrium 
radiation. Due to the differences in the exterior geome- 
try, the thermodynamic and electromagnetic properties 
of the disks (energy flux, temperature distribution and 
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equilibrium radiation spectrum) are different for these 
two classes of compact objects, thus giving clear obser- 
vational signatures, which may allow to discriminate, at 
least in principle, gravastars from black holes. We would 
like to point out that the proposed method for the detec- 
tion of the gravastars by studying accretion disks is an 
indirect method, which must be complemented by direct 
methods of observation of the surface of the considered 
compact objects. 

The present paper is organized as follows. In Sec. QTJ 
we present the fundamental field equations for static and 
slowly rotating gravastar models. In Sec. IIII1 we review 
the formalism and the physical properties of the thin 
disk accretion onto compact objects, for stationary ax- 
isymmetric spacetimes. In Sec. IIV1 we analyze the basic 
properties of matter forming a thin accretion disk around 
slowly rotating gravastar spacetimes. We discuss and 
conclude our results in Sec. [V] Throughout this work, 
we use a system of units so that c — G = K = kg = 1, 
where kg is Boltzmann's constant. 



II. SLOWLY ROTATING GRAVASTAR AND 
KERR BLACK HOLES 

In order to construct slowly rotating gravastar mod- 
els we first consider the static case. Then, by assuming 
that rotation represents a second order perturbation of 
the static slowly rotating gravastar model can be 

constructed. 



A. Static gravastar models 

1. Equations of structure 

For a static general relativistic spherically symmetric 
matter configuration, the interior line element can be 
taken generally as 

ds 2 = -e^dt 2 + e^dr 2 + r 2 {d9 2 + sin 2 Odq?) . (1) 

We assume that the star consists of an anisotropic fluid 
distribution of matter, and is given by 



(p + p±_)U^ U v + p± + (p r - p±_)x^Xv , (2) 



where C/' 1 is the four-velocity, is the unit spacelike 



vector in the radial direction, i.e., x p 



-A(r)/2 p 



p(r) is the energy density, p r (r) is the radial pressure 
measured in the direction of x^i and p±(r) is the trans- 
verse pressure measured in the orthogonal direction to 
Taking into account the above considerations, the 
stress-energy tensor is given by the following profile: 
T^v = dia,g[- p(r),p r {r),p ± (r),p ± (r)}. 

We suppose that inside the star p r ^ p±, Vr ^ 0, and 
define the anisotropy parameter as A = p± — p r , where 
A is a measure of the deviations from isotropy. If A > 
0, Vr the body is tangential pressure dominated while 



A < indicates that p r > p±. Note that A/r represents 
a force due to the anisotropic nature of the stellar model, 
which is repulsive, i.e., being outward directed if p± > p r , 
and attractive if p± < p r . 

The properties of the anisotropic compact object can 
be completely described by the gravitational structure 
equations, which are given by: 



dp r 
dr 



dm 2 

= iTrpr 

dr 



(p + p r ) \m + 4:TTr 3 p r 



2A 



dv 
dr 



(1-*?) 



dp r 



4A 



p + p r dr r (p + p r ) ' 



(3) 



(4) 



(5) 



where m(r) is the mass inside radius r, and the relation- 
ship e~ A ( r - ) = [1 — 2m(r)/r] has been used. 

A solution of Eqs. ©-© is possible only when bound- 
ary conditions have been imposed. As in the isotropic 
case we require that the interior of any matter distribu- 
tion be free of singularities, which imposes the condition 
m(r) — > as r — ► 0. Assuming that p r is finite at r = 0, 



we have 



as r — * 0. Therefore the gradient dp r /dr 



will be finite at r = only if A vanishes at least as rapidly 
as r when r — > 0. This requires that the anisotropy pa- 
rameter satisfies the boundary condition 



lim 



A(r) 



= 0. 



(6) 



At the center of the star the other boundary conditions 
for Eqs. ©-© are p r (0) = pj_(0) = p c and p(0) = p c , 
where p c and p c are the central density and pressure, 
respectively. The radius a of the star is determined by 
the boundary condition p r (a) = 0. We do not neces- 
sarily require that the tangential pressure p± vanishes 
for r — a. Therefore at the surface of the star the 
anisotropy parameter satisfies the boundary condition 
A (a) = p±{a) — p r (a) = p±(a) > 0. To close the field 
equations the equations of state of the radial pressure 
p r = p r (p) and of the tangential pressure p± = p^ (p) 
must also be given. To be a gravastar model, we need to 
impose the equation of state p r — —p, so that from the 
field equations one easily deduces that 



v(r) — — A(r) = In 



1 



2m(r) 



(7) 



where m(r) is the mass function. 



2. Junction interface 

We consider models of gravastars by matching an inte- 
rior solution, governed by an equation of state, p r = —p, 
to an exterior Schwarzschild vacuum solution with p = 
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p = 0, at a junction interface S, with junction radius a. 
The Schwarzschild metric is given by 



ds 2 



2M 



dt- 



r J 



dr z + r 2 dVL 2 



(8) 

which possesses an event horizon at r& = 2M, and dtt 2 = 
dO 2 + sin 2 9 dip 2 . To avoid the event horizon, the junction 
radius lies outside 2M, i.e., a > 2M. We show below 
that M, in this context, may be interpreted as the total 
mass of the gravastar. 

Using the Darmois-Israel formalism |34l|. the surface 
stresses of the thin shell are given by [1JJ 



I I 1 — — +d 2 +ad 1 — 771 ' — — +d 2 +aa 



P = 



(10) 



where the overdot denotes a derivative with respect to 
proper time, r. a and V are the surface energy density 
and the tangential pressure, respectively [ll|- The dy- 
namical stability of the transition layer of these compact 
spheres to linearized spherically symmetric radial pertur- 
bations about static equilibrium solutions was explored 
using Eqs. (f9|)- (fT0|) (see pjj for details). Large stability 
regions were found that exist sufficiently close to where 
the event horizon is expected to form, so that it would 
be difficult to distinguish the exterior geometry of these 
gravastars from astrophysical black holes. 

The surface mass of the thin shell is given by m s = 
Aim 2 a. By rearranging Eq. ©, evaluated at a static 
solution do, i.e., a = d = 0, one obtains the total mass of 
the gravastar, given by 



M = m(ao) + m s (ao) 



'1- 



2m(a ) m s (a ) 



2a n 



(11) 



where bo is a non- negative constant, which was exten- 
sively analyzed in [ill ] in the context of dark energy stars. 
The latter may be determined from the regularity con- 
ditions and the finite character of the energy density at 
the origin r — 0, and is given by &o = 8irp c /3, where p c 
is the energy density at r = 0. 

This choice of the mass function represents a mono- 
tonic decreasing energy density in the star interior, and 
was used previously in the analysis of isotropic fluid 
spheres by Matese and Whitman |35l| as a specific case 
of the Tolman type— TV solution [36}, and later by Finch 
and Skea [37]] • Anisotropic stellar models, with the re- 
spective astrophysical applications, were also extensively 
analyzed in Refs. (38|, by considering a specific case of 
the Matese- Whitman mass function. The numerical re- 
sults outlined show that the basic physical parameters, 
such as the mass and radius, of the model can describe 
realistic astrophysical objects such as neutron stars [38|. 

The spacetime metric for this solution is provided by 



ds 2 = - 



{^)dt 2 +(^)dr 2 + Pdn 2 .(U) 



The stress-energy tensor components are given by 



Pr 



b \ (3 + 26 r 2 ) 



P 



\8ttJ {l + 2b r 2 ) 2 
b \ (3 + 2b r 2 ) fb 2 r 2 \ (5 + 2b r 2 ) 



8tt J (1 + 2b Q r 2 ) 2 V 47r / (1 + 2fe r 2 ) 3 
The anisotropy factor takes the following form 
b 2 r 2 \ (5 + 2& r 2 ) 



A = 



4tt 7(1 + 2b r 2 f 



(14) 



which is always positive, implying that p± > p r , and 
A| r=0 = at the center, i.e., p±(0) = p r (0), as expected. 



3. Specific model: Tolman- Matese- Whitman mass function 



B. Slowly rotating gravastars 



To gain insight into the problem, it is interesting to 
present a specific example. For instance, consider the 
following choice for the mass function, given by 



m(r) 



b r 3 



2(l + 2b r 2 ) 



(12) 



When the equilibrium configuration described in Sec. 
Ill A H is set into slow rotation, the geometry of spacetime 
around it and its interior distribution of stress-energy are 
changed. With an appropriate change of coordinates, the 
perturbed geometry is described by [1^, |40| ] 



ds 2 



{1 + 2[h + h 2 P 2 (cos 9)}} dt 2 



1 + 2 [rn + m 2 P2 (cos 9)} I (r - 2M) 



1 - 2M/r 

[1 + 2 (v 2 -h 2 )P 2 (cos 6>)] \d9 2 + sin 2 9 (dip - tudt) 2 ] + O (fl%) 



dr 2 



(15) 



Here, P 2 (cos 6*) = (3 cos 2 9 — l) /2 is the Legendre polynomial of order two; the angular velocity, w, of the local 
inertial frame, is a function of the radial coordinate r, and is proportional to the star's angular velocity fig; and ho, 
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/i2, mo, 771-2 , V2 are functions of r that are proportional to Oi. Outside the star the metric can be written as [4C 



ds 2 



2M J 2 
1 + 2— 

r 

2M J 2 
1 + 2— 

2 

Mr 3 



1 + 2 



1 - 2 



J 

Mr 3 

' J 2 



r 



5 Q-jyM 2 
8 M3 ^ 2 



(X) 



Mr 3 



5M 



J 2 /M 



j 1 + 2 

x J d6» 2 + sin 2 



8 M 3 ^ 2 (X) 
2M 



P 2 (cos 6>) | dt 2 
P 2 (cos6»)lrfr 2 



J 2 / 2M\ 5Q-J 2 /M . . . 

tt q ( i + — I + - , r3 ' ( ; = — Ql (x) - Ql (x) 



2J 



r / 

2 ' 



M 3 



y/r (r - 2M) 



P 2 (cost 



dip jdt 



(16) 



In Eq. (p~6j) the variable x = r/M—1, and the quantities 
Q2 an d <?! denote associated Legendre polynomials of the 
second time, so that 



Ux) = Vx^T 



(3x 2 



and 



(X) 



5x ~ 3x 2 

(x 2 - 1) 



(x 2 



3 . X+l 

oXln 7 

2 X- 1 



1) In- 



x-i' 



(17) 



(18) 



respectively. The line element outside the star is deter- 
mined by three constants: the total mass of the rotating 
star M, the star's total angular momentum J and the 
star's mass quadrupole moment Q. The quadrupole mo- 
ment Q can be expressed in terms of the eccentricity e = 

\J(r e /r p ) 2 - 1 of the star as e = y/3Q/Mr* 2 + (l/r* 2 ), 
where r* is a large distance from the origin, and r e and 
r p are the equatorial and polar radius of the star, re- 
spectively [40]. The metric given by Eq. (fT6|) is valid 
in the case of slow rotation, that is, the angular veloc- 
ity of the star fig must be small enough so that the 
fractional changes in pressure, energy density, and grav- 
itational field, due to rotation, are all much less than 
unity. The condition of slow rotation can be formulated 
as Q| < (c/Rf (GM/Rc 2 ), where R is the radius of 
the static stellar configuration [3!|. The critical angu- 
lar velocity at which mass shedding occurs is given by 
Q K = GM/R 3 [lO], and therefore the condition of slow 
rotation implies fls <C VIk ■ All models considered in the 
present paper satisfy this condition. 

The metric given by Eq. (|16[) is used to determine 
the electromagnetic signatures of accretion disks around 
slowly rotating gravastars, which is analyzed in detail 
below. 



The Kerr metric, describing a rotating black hole, is given 
in the Boyer-Lyndquist coordinate system by 



ds z 



1 



2Mr\ 



dt 



2Mr 



-a sin 8dtd(f> 



^ K -dr 2 



2Mr 



+Y. K d9 2 + yr 2 + a 2 + a 2 sin 2 9j sin 2 9d<p 2 , (19) 

where = r 2 + a 2 cos 2 9 and Ak = r 2 + a 2 — 2mr, 
respectively. In the equatorial plane, the metric compo- 
nents reduce to 



9tt = 



9t<p 



Qrr 

9<f><P 



-(' 


2Mr 


£k 


2Mr 






a sin 2 9 




r 2 


A K 


" A K ' 




2 2 
-a 2 + - 



.Ma 



2 ■ 2 
-a sin 1 



sin 



Sk 
2M\ 
—) ' 

respectively. For the Kerr metric J = —Ala and 
Q = J 2 /M, respectively. The latter relationship, i.e., 
Q = J 2 /M , between the quadrupole momentum and the 
angular momentum shows the very special nature of the 
Kerr solution. 



III. ELECTROMAGNETIC RADIATION 
PROPERTIES OF THIN ACCRETION DISKS IN 
STATIONARY AXISYMMETRIC SPACETIMES 

To set the stage, we present the general formalism 
of electromagnetic radiation properties of thin accretion 
disks in stationary axisymmetric spacetimes. 



C. The Kerr black hole 

For self-completeness and self-consistency, we present 
the Kerr metric, as it will be compared to metric (|16|) in 
the electromagnetic signature analysis of accretion disks. 



A. Stationary and axially symmetric spacetimes 

In this work we analyze the physical properties and 
characteristics of particles moving in circular orbits 
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around general relativistic compact spheres in a station- 
ary and axially symmetric geometry given by the follow- 
ing metric 



. dr 2 + 



,de 2 



ds 2 = gtt dt 2 + 2g t< f, dtd<j> + g r 

(20) 

Note that the metric functions gtt, gt<f>, 9rr, gee and 
only depend on the radial coordinate r in the equatorial 
approximation, i.e., \6 — ir\ <C 1, which is the case of 
interest here. In the following we denote the square root 
of the determinant of the metric tensor by \J— g. 

To compute the flux integral given by Eq. (|39[) , we de- 
termine the radial dependence of the angular velocity Q, 
of the specific energy E and of the specific angular mo- 
mentum L of particles moving in circular orbits around 
compact spheres in a stationary and axially symmetric 
geometry through the geodesic equations. The latter take 
the following form 



9r 



dt 

d<p 

2 



Eg </></> + Lg t<j> 
9% ~ 9tt9<t><t> ' 

Egtd, + Lgtt 



- 9tt9W 

E 2 g u + 2ELg tl)) + L 2 g tt 
9^ 



9ug4,<p 

One may define an effective potential term defined as 
E 2 g<h<h + 2ELg trt , + L 2 g tt 



(21) 
(22) 
(23) 



V ef f(r) = -1 + 



9t<t> 



9ugu 



(24) 



For stable circular orbits in the equatorial plane 
the following conditions must hold: V e ff(r) = and 
Veff, r(f) = 0. These conditions provide the specific en- 
ergy, the specific angular momentum and the angular 
velocity of particles moving in circular orbits for the case 
of spinning general relativistic compact spheres, given by 



E = 
L = 

n = 



gtt + gt<t>Q 



\/-gtt - 2^0 

9t<t> - 



' 94>4> 



n 2 



\f-9tt ~ 2gt<t>n 

d(j) _ —9t(j>,r + \/(gt0,r) 2 ~ gtt,rg00,r 
dt Qchch.r 



(25) 
(26) 

(27) 



The marginally stable orbit around the central object can 
be determined from the condition V r e // i rr (r) = 0, which 
provides the following important relationship 



E 2 g^,; 



+ 2ELgt<j } , rr + L 2 g t t ir r — (5t0 — 9ttg<j><j>) ,rr — 0- 



(28) 



By inserting Eqs. (gSJ-lgTl) into Eq. fl55]) and solving this 
equation for r, we obtain the marginally stable orbit for 
the explicitly given metric coefficients gtt, gt<p and g^. 



For a Kerr black hole the geodesic equation ([23]) for r 
becomes 



(29) 



with the effective potential given by 



V eff (r) = -1 + \E 2 [r 2 {r 2 + a 2 ) + 2ma 2 



-AELmar - L 2 



2mr 



}/ V 2 



2mr 



2 )} 



(30) 



Note that these relationships may be rewritten in the 
following manner 



1Z) =V(r) 



(31) 



with V(r) given by 

V(r) =r 2 A K V ef f(r) = r 2 (r 2 - 2mr + a 2 )V e} f{r) . (32) 

where the relationship Ak — g 2 ^ — 9ttg<p4> = r 2 — 2mr + a 2 
along the equatorial plane has been used. 

B. Properties of thin accretion disks 

For the thin accretion disk, it is assumed that its verti- 
cal size is negligible, as compared to its horizontal exten- 
sion, i.e, the disk height H , defined by the maximum half 
thickness of the disk in the vertical direction, is always 
much smaller than the characteristic radius r of the disk, 
defined along the horizontal direction, H <C r. The thin 
disk is in hydrodynamical equilibrium, and the pressure 
gradient and a vertical entropy gradient in the accret- 
ing matter are negligible. The efficient cooling via the 
radiation over the disk surface prevents the disk from cu- 
mulating the heat generated by stresses and dynamical 
friction. In turn, this equilibrium causes the disk to sta- 
bilize its thin vertical size. The thin disk has an inner 
edge at the marginally stable orbit of the compact ob- 
ject potential, and the accreting plasma has a Keplerian 
motion in higher orbits. 

In steady state accretion disk models, the mass ac- 
cretion rate Mo is assumed to be a constant that does 
not change with time. The physical quantities describ- 
ing the orbiting plasma are averaged over a characteristic 
time scale, e.g. At, over the azimuthal angle A(f> = 2ir 
for a total period of the orbits, and over the height H 

The particles moving in Keplerian orbits around the 
compact object with a rotational velocity Q — d<p/dt have 
a specific energy E and a specific angular momentum L, 
which, in the steady state thin disk model, depend only 
on the radii of the orbits. These particles, orbiting with 
the four- velocity form a disk of an averaged surface 
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density S, the vertically integrated average of the rest 
mass density pg of the plasma. The accreting matter in 
the disk is modeled by an anisotropic fluid source, where 
the density po , the energy flow vector and the stress 
tensor t^ v are measured in the averaged rest-frame (the 
specific heat was neglected). Then, the disk structure 
can be characterized by the surface density of the disk 

MM, 



can be expressed in terms of the specific energy, angu- 
lar momentum and of the angular velocity of the central 
compact object [22I l24j] . 



F{r) 



M 



n 



my 



(E-QL)L, r dr . (39) 



E(r) 



H 



(Po)dz, 



(33) 



—H 



with averaged rest mass density (po) over At and 2n and 
the torque 



H 



Wi = / (U r )dz, 



(34) 



H 



with the averaged component (t£) over At and 2ir. The 
time and orbital average of the energy flux vector gives 
the radiation flux T(r) over the disk surface as T{r) — 
{q z )- 

The stress-energy tensor is decomposed according to 



= p^u" + 2u {fl q u) + f", 



(35) 



where = 0, u^t^ = 0. The four-vectors of 

the energy and angular momentum flux are defined by 
-E^ = Tl l {d/dty and J' 1 = T£{d/d4>) v , respectively. 
The structure equations of the thin disk can be derived 
by integrating the conservation laws of the rest mass, of 
the energy, and of the angular momentum of the plasma, 
respectively [12, HH . From the equation of the rest mass 
conservation, V AI (pow' i ) = 0, it follows that the time aver- 
aged rate of the accretion of the rest mass is independent 
of the disk radius, 



Mn 



-2 r K^/—gYiU r = constant. 



(36) 



Another important characteristics of the mass accre- 
tion process is the efficiency with which the central object 
converts rest mass into outgoing radiation. This quan- 
tity is defined as the ratio of the rate of the radiation of 
energy of photons escaping from the disk surface to infin- 
ity, and the rate at which mass-energy is transported to 
the central compact gene ral relativistic object, both mea- 
sured at infinity [UHil- If all the emitted photons can 
escape to infinity, the efficiency is given in terms of the 
specific energy measured at the marginally stable orbit 



l-E n 



(40) 



For Schwarzschild black holes the efficiency e is about 
6%, whether the photon capture by the black hole is con- 
sidered, or not. Ignoring the capture of radiation by the 
hole, e is found to be 42% for rapidly rotating black holes, 
whereas the efficiency is 40% with photon capture in the 
Kerr potential [25[. 

The accreting matter in the steady-state thin disk 
model is supposed to be in thermodynamical equilibrium. 
Therefore the radiation emitted by the disk surface can 
be considered as a perfect black body radiation, where 
the energy flux is given by F(r) = <TsBT 4 (r) (ctsb is the 
Stefan-Boltzmann constant), and the observed luminos- 
ity L (y) has a redshifted black body spectrum [28[ : 



The conservation law V^-E^ = of the energy has the 
integral form 



[M E - 2iry/=gSlW/], 



-gFE 



(37) 



which states that the energy transported by the rest mass 
flow, MqE, and the energy transported by the dynami- 
cal stresses in the disk, 2-k^/— gflW^ , is in balance with 
the energy radiated away from the surface of the disk, 
An^—gFE. The law of the angular momentum conser- 
vation, V M J M = 0, also states the balance of these three 
forms of the angular momentum transport, 



[M L - 27rrWy] ir = 4iry/=gFL 



(38) 



By eliminating W<p r from Eqs. (|37j) and (|38)l . and 
applying the universal energy-angular momentum rela- 
tion dE = fldJ for circular geodesic orbits in the form 
E tr = ilL yr , the flux F of the radiant energy over the disk 



L (v) = And I (y) = - cos 7 

7T 



exp [v e /T) - I 



(41) 

Here d is the distance to the source, I{v) is the Planck 
distribution function, 7 is the disk inclination angle, and 
Ti and rj indicate the position of the inner and outer 
edge of the disk, respectively. We take r,j = r ms and 
rj — > 00, since we expect the flux over the disk surface 
vanishes at r — > 00 for any kind of general relativistic 
compact object geometry. The emitted frequency is given 
by v e = v(l + z), where the redshift factor can be written 
as 



1 + z 



1 + fir sin 4> sin 7 



\J-9tt - 2ilg 



t<t> 



(42) 



where we have neglected the light bending [4l|, |42 
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IV. ELECTROMAGNETIC AND 
THERMODYNAMIC SIGNATURES OF 
ACCRETION DISKS AROUND SLOWLY 
ROTATING GRAVASTARS AND BLACK HOLES 

A. Electromagnetic and thermodynamic properties 
of the disks 

In this section we compare the radiation properties of 
thin accretion disks around gravastars and black holes in 
the slowly rotating case when the spin parameter a* — 
J/M 2 has the maximal value of 0.5. In Fig. [1] we present 
the time averaged energy flux F(r) radiated by the disk 
for both types of central objects with the total mass M 
of 10 6 M Q and increasing spin parameter from 0.1 to 0.5. 
The quadrupole moment Q of the gravastar models runs 
between 0.1M 3 and 2M 3 . Here the mass accretion rate 
M is set to 2.5 x 10 _5 A/©/yr, which is in the typical 
range for super massive cental objects. 

If we compare the flux emerging from the surface of the 
thin accretion disk around black holes and gravastars, 
we find that its maximal value is systematically lower 
for gravastars, independently of the values of the spin 
parameter or the quadrupole momentum. For very slow 
rotation, (a* = 0.1), and a relative small value of the 
quadrupole moment (Q = 0.1M 3 ), the radial distribution 
of the disk radiation is close to each other for the two 
types of compact central objects. The maximal flux for 
gravastars is roughly 90% of the black hole's flux, and 
the maximum of the inner edge of the accretion disk is 
located at somewhat higher radii for gravastars. With 
increasing rotational frequency of the central object, the 
flux values also increase, but the increment is higher for 
black holes than for gravastars. For a* = 0.5 the flux 
maximum for black holes is almost twice the maximal flux 
value for gravastars. The more rapid rotation does not 
cause strong effect on the location of the inner disk edge 
for gravastars, as we find a slight decrease in the value 
of r ms as the rotation of central object increases (see in 
Tab HI)). For black holes, the effect of the rotation is also 
stronger here. Although the error of the approximation 
applied for the slow rotation is rapidly increasing in the 
regime around a* ~ 0.5, this picture is still definitely 
adequate for lower values of a*. 

The variation of the quadrupole moment causes con- 
siderable changes in both the maximal value of disk ra- 
diation and the location of the inner edge of the disk. 
As we increase Q, the maximal flux decreases and r ms 
increases. These effects are presented in Fig. [21 show- 
ing the disk temperatures, although the differences are 
somewhat less striking. 

The disk spectra, presented in Fig. [3J have similar fea- 
tures in the dependence of the disk radiation on the ro- 
tation parameter and the quadrupole momentum. The 
amplitudes and the cut-off frequencies of the spectra for 
gravastars are always lower than those for black holes. 
For higher rotational velocity, the amplitudes are some- 
what higher but do not exhibit much change. The cut-off 



frequency for black holes increases moderately, whereas 
it has only a negligible increment for gravastars. This 
makes the differences in the spectral properties more 
acute for higher values of the spin parameter (a* ^ 0.3). 
The increase in the quadrupole moment somewhat low- 
ers the amplitude of the spectra but causes a stronger 
decrease in the cut-off frequencies. 



B. Conversion efficiency of the accreting mass 

We also present the conversion efficiency e of the ac- 
creting mass into radiation, measured at infinity, which 
is given by Eq. ([40)) , for the case where the photon cap- 
ture by the slowly rotating central object is ignored. The 
value of e measures the efficiency of energy generating 
mechanism by mass accretion. The amount of energy 
released by matter leaving the marginally stable orbit, 
and falling down the black hole, is the binding energy 
E ms of the black hole potential. In Tabs. HI and HP the 
marginally stable orbits r ms and e are given for black 
holes and gravastars with the parameters a* and Q in the 
range used for the plots presenting the radiation proper- 
ties of the accretion disks. 





Vms [M] 


t [io- 2 ] 


0.1 


5.67 


6.06 


0.2 


5.33 


6.46 


0.3 


4.98 


6.94 


0.4 


4.62 


7.51 


0.5 


4.24 


8.21 



TABLE I: The marginally stable orbit and the efficiency for 
slowly rotating Kerr black holes with different spin parame- 
ters. 

These values demonstrate the variation in the location 
of the inner disk edge with the changing spin parame- 
ter and quadrupole momentum, as we have seen in the 
discussion on the radial distribution of the flux. The 
higher these values are, the closer the marginally sta- 
ble orbits are to the center in the dimensionless radial 
scale. For slowly rotating black holes, the conversion ef- 
ficiency is still close to 6%, which is the value obtained 
for Schwarzschild black holes. Up to a* = 0.5 it increases 
to 8%, which is still much lower than the one for extreme 
Kerr black holes. For gravastars, e has a smaller vari- 
ation, and always remains smaller than the conversion 
efficiency for black holes. For very slow rotation, the ef- 
ficiency is approximately 5.8%, which is close to the one 
for the static black hole, and it decreases to 5% as the 
quadrupole moment increases to 2M 3 . For a higher spin 
parameter (a* ~ 0.5), e is still around 6.5% but it be- 
comes smaller than 6% as Q increases. We conclude that 
the conversation efficiency is higher for more rapidly ro- 
tating gravastars, but this is moderated by the increment 
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FIG. 1: The energy flux emerging from the accretion disk of slowly rotating gravastars and black holes for the spin parameter 
a* = 0.1 (upper left hand plot), a* = 0.3 (upper right hand plot), a, = 0.4 (lower left hand plot), and a, = 0.5 for (lower right 
hand plot). All the plots are given for the total mass M = 1O 6 M0, the quadrupole moments Q — 0.1,0.5, 1.0,2.0 times M 3 , 
and the mass accretion rate 2.5 x 10~ 5 M©/yr. 




FIG. 2: The disk temperature of slowly rotating gravastars and black holes for the spin parameter a* = 0.1 (upper left hand 
plot), a» = 0.3 (upper right hand plot), a» = 0.4 (lower left hand plot), and a* = 0.5 for (lower right hand plot). All the plots 
are given for the total mass M = 1O 6 M0, the quadrupole moments Q = 0.1, 0.5, 1.0, 2.0 times M 3 , and the mass accretion rate 
2.5 x lO~ 5 M /yr. 
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v [Hz] v [Hz] 



FIG. 3: The disk spectra of slowly rotating gravastars and black holes for the spin parameter a* =0.1 (upper left hand plot), 
a, — 0.3 (upper right hand plot), a* = 0.4 (lower left hand plot), and a* = 0.5 for (lower right hand plot). All the plots are 
given for the total mass M = 1O 6 M0, the quadrupole moments Q — 0.1,0.5,1.0,2.0 times M 3 , and the mass accretion rate 
2.5 x 10~ 5 M Q /yr. 



a* 


Q [M 3 ] 


r m3 [M] 


e [10- 2 ] 


0.1 


0.1 


5.91 


5.82 


0.1 


0.5 


6.20 


5.60 


0.1 


1.0 


6.50 


5.37 


0.1 


2.0 


7.00 


5.02 


0.2 


0.1 


5.75 


6.00 


0.2 


0.5 


6.05 


5.75 


0.2 


1.0 


6.37 


5.50 


0.2 


2.0 


6.77 


5.23 


0.3 


0.1 


5.58 


6.19 


0.3 


0.5 


5.89 


5.90 


0.3 


1.0 


6.23 


5.63 


0.3 


2.0 


6.89 


5.12 


0.4 


0.1 


5.39 


6.40 


0.4 


0.5 


5.74 


6.08 


0.4 


1.0 


6.09 


5.77 


0.4 


2.0 


6.65 


5.33 


0.5 


0.1 


5.20 


6.64 


0.5 


0.5 


5.58 


6.27 


0.5 


1.0 


5.95 


5.93 


0.5 


2.0 


6.52 


5.45 



TABLE II: The marginally stable orbit and the efficiency for 
slowly rotating gravastars with different spin parameters and 
quadrupole moments. 



in its quadrupole moment. In addition to this, it is al- 
ways smaller than the e for black holes, i.e., gravastars 
provide a less efficient mechanism for converting mass to 
radiation than black holes. 

In order that our proposal for discriminating gravastars 
from black holes by using the electromagnetic emission 
properties of accretion disks could be effectively applied 
in concrete observational cases, it is necessary to know 
at least the value of the mass and of the spin parameter 
of the central rotating compact object. To estimate the 
spins of stellar-mass black holes in X-ray binaries, one 
has to fit the continuum X-ray spectrum of the radia- 
tion from the accretion disk, by using the standard thin 
disk model, and extract the dimensionless spin parame- 
ter a* = a/M of the black hole as a parameter of the fit 
[44| . Recently, a number of precise black hole spin de- 
terminations have been reported. By using Chandra and 
Gemini-North observations of the eclipsing X-ray binary 
M33 X-7, precise values of the mass of its black hole pri- 
mary and of the system's orbital inclination angle have 
been obtained. The distance to the binary is also known 
to a few percent. By using these precise results, from the 
analysis of 15 Chandra and XMM-Newton X-ray spectra, 
and a fully relativistic accretion disk model, one can find 
that the dimensionless spin parameter of the black hole 
primary is a* = 0.77 ± 0.05 [45[. Therefore, even that 
presently there are severe observational limitations for 
the application of our proposal, with the future improve- 
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merits of the observational techniques, the observation 
of the emission spectra of accretion disks could be effec- 
tively used to discriminate between gravastars and black 
holes. 



V. DISCUSSIONS AND FINAL REMARKS 

If gravastars are surrounded by a thin shell of mat- 
ter, the presence of a turning point for matter (the point 
where the motion of the infalling matter suddenly stops) 
at the surface of the gravastars may have important as- 
trophysical and observational implications. Since the 
velocity of the matter at the gravastar surface is zero, 
matter can be captured and deposited on the surface of 
gravastar. Therefore gravastars may have a gaseous sur- 
face, formed from a a thin layer of dense and hot gas. 
Moreover, because matter is accreted continuously, the 
increase in the size and densit y o f the surface will ignite 
some thermonuclear reactions [29[ . The ignited reactions 
are usually unstable, causing the accreted layer of gas 
to burn explosively within a very short period of time. 
After the nuclear fuel is consumed, the gravastar reverts 
to its accretion phase, until the next thermonuclear in- 
stability is triggered. Thus, gravastars may undergo a 
semi-regular series of explosions, called type I thermonu- 
clear bursts, discovered first for X-ray binaries (46l. |47|. 

The observational signatures indicating the presence 
of X-ray bursts from gravastars are similar to those of 
standard neutron stars, and are the gravitational redshift 
of a surface atomic line, the touchdown luminosity of a 
radius-expansion burst, and the apparent surface area 
during the cooling phases of the burst [48[ . 

If the thermal radiation with wavelength A e emitted by 
the matter at the surface of the gravastar has absorption 
or emission features characteristic of atomic transitions, 
these features will be detected at infinity at a wavelength 
A D , gravitationally redshifted with a value 



A — A e 

Ae 



= e -Kfl)/2 _ L) 



(43) 



where we have assumed, for simplicity, that the gravas- 
tar is static, and that the exterior metric can be approx- 
imated by the standard Schwarzschild metric, given by 
Eq. $5§. By assuming a gravastar of mass M = 4x 10 6 M Q 
and radius B = 1.4 x 10 12 cm, we obtain a surface red- 
shift of z grav — 1.55. The corresponding value of the 
redshift for a neutron star with mass M = 2M and ra- 
dius R = 10 6 cm is zns = 0.56. Therefore the radiation 
coming from the surface of a gravastar may be highly red- 
shifted (in standard general relativity the redshift obeys 
the constrain z < 2). 

Type I X -ray bursts show strong spectroscopic evi- 
dence for a rapid expansion of the radius of the X-ray 
photosphere. The luminosities of these bursts reach the 
Eddington critical luminosity at which the outward radi- 
ation force balances gravity, causing the expansion lay- 
ers of the star. The touchdown luminosity of radius- 



expansion bursts from a given source remain constant 
between bursts to within a few percent, giving empirical 
verification to the theoretical expectation that the emerg- 
ing luminosity is approximately equal to the Eddington 
critical luminosity. The Eddington luminosity at infinity 
of a gravastar is given by [48| 



47rm i? 2 _ A 



de»' 2 



dr 



(44) 



r=R 



where mo is the mass of the particle and a is the inter- 
action cross section. For the gravastar we obtain 



roo 

L E — 



4irm M 



2M 
r 



(45) 



The ratio of the Eddington luminosities at infinity for 
a gravastar with mass of 4 x 10 6 Mq and radius 1.4 x 10 12 
cm and a neutron star with mass of 2 solar masses and 
radius of 10 km is Lf grav /Lf NS = 1.22 x 10 6 . Finally, 
we consider the apparent surface area during burst cool- 
ing. Observations of the cooling tails of multiple type 
I bursts from a single source have shown that the ap- 
parent surface area of the emitting region, defined as 
S°° = 4 / kD 2 F CjO0 /(jsbT^< x> , where F C)00 is the measured 
flux of the source during the cooling tail of the burst, 
T Ci oo is the measured color temperature of the burst spec- 
trum, D is the distance to the source and asB is the 
Stefan-Boltzmann constant, remains approximately con- 
stant during each burst, and between bursts from the 
same source. The color temperature on the surface of 
the compact object T c ^ is related to the color temper- 



ature measured at infinity by T c . 



By introducing the color correction factor f c — T c /T e ff, 
where T e f f is the effective temperature at the surface, we 
obtain 



S°° = 4^[z(R) + lf 

J C 



(46) 



Since the radius of the gravastar as well as its surface 
redshift may be very large quantities, the apparent area 
of the emitting region as measured at infinity may be also 
very large. Hence all the astrophysical quantities related 
to the observable properties of the X-ray bursts, origi- 
nating at the surface of the gravastar can be calculated, 
and have finite values on the surface of the gravastar and 
at infinity. 

It was argued in [291 ] that any neutron star, composed 
by matter described by a more or less general equation 
of state, should experience thermonuclear type I bursts 
at appropriate mass accretion rates. The question asked 
in [29j is whether an "abnormal" surface may allow such 
a behavior. The gravastars may have such a zero veloc- 
ity, particle trapping, abnormal surface. The presence 
of a material surface located at the "event horizon" im- 
plies that energy can be radiated, once matter collides 
with that surface. Thus, gravastar models, characterized 
by high mass, normal matter crusts/surfaces and type 
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I thermonuclear bursts can be theoretically constructed. 
Moreover, some of the so-called SXT's (soft X-ray tran- 
sients), having a relatively low mass function ( e.g . SXT 
A0620-00, with a mass function f(M) > 3M Q HI), but 
still exceeding the equilibrium limit of 3M Q , or very mas- 
sive neutron stars showing the presence of a crust, may 
in fact be gravastars. 

It is generally expected that most of the astrophysical 
objects grow substantially in mass via accretion. Re- 
cent observations suggest that around most of the active 
galactic nuclei (AGN's) or black hole candidates there 
exist gas clouds surrounding the central far object, and 
an associated accretion disk, on a variety of scales from 
a tenth of a parsec to a few hundred parsecs [43| • These 
clouds are assumed to form a geometrically and optically 
thick torus (or warped disk), which absorbs most of the 
ultraviolet radiation and the soft x-rays. The gas exists in 
either the molecular or the atomic phase. Evidence for 
the existence of super massive black holes comes from 
the very long baseline interferometry (VLBI) imaging of 
molecular H 2 masers in active galaxies, like NGC 4258 
and from the astrometric and radial velocity mea- 
surements of the fully unconstrained Keplerian orbits for 
short period stars around the supermassive black hole at 
the center of our galaxy \EM- The VLBI imaging, 
produced by Doppler shift measurements assuming Kep- 
lerian motion of the masering source, has allowed a quite 
accurate estimation of the central mass, which has been 
found to be a 3.6 x 10 7 M© super massive dark object, 
within 0.13 parsecs. Hence, important astrophysical in- 
formation can be obtained from the observation of the 
motion of the gas streams in the gravitational field of 
compact objects. 

Therefore the study of the accretion processes by com- 
pact objects is a powerful indicator of their physical na- 
ture. However, up to now, the observational results have 
confirmed the predictions of general relativity mainly in 
a qualitative way. With the present observational preci- 
sion one cannot distinguish between the different classes 
of compact/exotic objects that appear in the theoreti- 
cal framework of general relativity [29(. However, with 
important technological developments one may allow to 
image black holes and other compact objects directly [52| . 
Recent observations at a wavelength of 1.3 mm have set a 
size of microarcseconds on the intrinsic diameter of Sgr A* 
[HI . This is less than the expected apparent size of the 
event horizon of the presumed black hole, thus suggest- 
ing that the bulk of SgrA* emission may not be centered 
on the black hole, but arises in the surrounding accre- 
tion flow. A model in which Sgr A* is a compact object 
with a thermally emitting surface was considered in [54j | . 
Given the very low quiescent luminosity of Sgr A* in 
the near-infrared, the existence of a hard surface, even 
in the limit in which the radius approaches the hori- 
zon, places a severe constraint on the steady mass ac- 
cretion rate onto the source: M < lO _12 M0/yr. This 
limit is well below the minimum accretion rate needed 
to power the observed submillimeter luminosity of Sgr 



A*: M > KT lo Af /yr. Thus it follows that Sgr A* 
does not have a surface, i.e., that it must have an event 
horizon. This argument could be made more restrictive 
by an order of magnitude with microarcsecond resolution 
imaging, e.g., with submillimeter very long baseline inter- 
ferometry. Submilliarcsecond astrometry and imaging of 
the black hole Sgr A* at the Galactic Centre may become 
possible in the near future at infrared and submillimetre 
wavelengths [55( ] . The expected images and light curves, 
including polarization, associated with a compact emis- 
sion region orbiting the central black hole were computed 
in ,56] . From spot images and light curves of the observed 
flux and polarization it is possible to extract the black 
hole mass and spin. At radio wavelengths, disc opacity 
produces significant departures from the infrared behav- 
ior, but there are still generic signatures of the black hole 
properties. Detailed comparison of these results with fu- 
ture data can be used to test general relativity, and to 
improve existing models for the accretion flow in the im- 
mediate vicinity of the black hole. 

With the improvement of the imaging observational 
techniques, it will also be possible to provide clear obser- 
vational evidence for the existence of gravastars, and to 
differentiate them from other types of compact general 
relativistic objects. Indeed, in this work we have shown 
that the thermodynamic and electromagnetic properties 
of the disks (energy flux, temperature distribution and 
equilibrium radiation spectrum) are different for these 
two classes of compact objects, consequently giving clear 
observational signatures. More specifically, comparing 
the energy flux emerging from the surface of the thin ac- 
cretion disk around black holes and gravastars of similar 
masses, it was found that its maximal value is system- 
atically lower for gravastars, independently of the val- 
ues of the spin parameter or the quadrupole momentum. 
These effects are confirmed from the analysis of the disk 
temperatures and disk spectra. In addition to this, it is 
also shown that the conversion efficiency of the accreting 
mass into radiation is always smaller than the conversion 
efficiency for black holes, i.e., gravastars provide a less ef- 
ficient mechanism for converting mass to radiation than 
black holes. Thus, these observational signatures may 
provide the possibility of clearly distinguishing rotating 
gravastars from Kerr- type black holes. 
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